Reply to Comment on "Critical analysis of a variational method used to describe 

molecular electron transport" 

loan BaldecQ and Horst Koppel 
Theoretische Chemie, Physikalisch-Chemisches Institut, Umversitat Heidelberg, 
Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany 
(Dated: August 5, 2011) 

We show that the failure of the Delaney-Greer (DG) variational ansatz for transport demonstrated 
by us in Phys. Rev. B 80, 165301 (2009) (I) is not related to an unsuitable constraint that prevents 
a broken time-reversal symmetry or to real orbitals, as DG incorrectly claim. The complex orbitals 
suggested by them as a way-out solution merely represent a particular case of the general case 
considered by us in I, which do not in the least affect our conclusion. In conjunction with the issues 
raised by the DG's Comment, we show that the DG Wigner conditions can erroneously constrain 
outgoing and not incoming charge carriers and present an example revealing that the sign of the 
"momentum" P of the Wigner function f[q, P) is not necessarily associated with the direction of 
motion in the real world. We also discuss a general reason why transport approaches which, like 
the DG's, are solely based on information on an isolated nanocluster are incorrect. 

I. INTRODUCTION 

In their Commentii Dclancy and Grccr (DG) claim that the unphysical result (current J = 0) presented by us in 
Ref.d (hereafter, I) as evidence against their variational approach^ is due to the fact (i) that our constraints (electronic 
populations) prevent a broken time-reversal symmetry, and (ii) that we used real orbitals, vifhich concomitantly 
constrain incoming and outgoing electrons. Further, DG accept that our idea to constrain populations "is an interesting 
alternative to the Wigner function" (WF), but (iii) argue that to get J 7^ it is essential to use complex orbitals, like 
those obtained by applying periodic boundary conditions (pBCs) for electrodes. By using the superscript DG and by 
referring to their work^ which we amply criticized^ii^i^ they aim to convey the false impression that they used such 
complex orbitals previously. 

Responding briefly, claim (i) is incorrect, claim (ii) has absolutely no real basis, while claim (iii) does not mender 
their variational approach; the complex orbitals discussed by them solely represent a particular case of the general 
case of I and do not in the least change the unphysical linear response current J = 0. 

II. REBUTTAL OF THE CRITICISM THAT POPULATION CONSTRAINTS PREVENT A BROKEN 

TIME-REVERSAL SYMMETRY 

DG argue that the unphysical result (J ~ 0) deduced in I is the consequence of our unsuitable choice of constraints, 
namely the particle number operators. They (incorrectly) claim that constraining real particle number operators does 
not break time-reversal symmetry, and for this very reason the current vanishes even beyond the linear response limit 
examined in I. 

We shall immediately show that this is not true: our constraints allow a broken time-reversal symmetry. The fact 
that this symmetry is not broken is the result of the defective DG variational ansatz (constrained energy minimization 
at zero temperature). Below, we shall only discuss the particular constraint — = a\aK, mentioned in the 
Comment. Let us consider a general many-body state, |^') = ^„ A„ |4'„), expanded in terms of the complete set of the 
real eigcnstates of the Hamiltonian H of the finite transport cluster without external bias {V — 0), -ff|5'„) = £'„|5'„). 
Notice that this exact expansion is just in the spirit of the DG approach, which is based on a configuration interaction 
(CI), albeit approximate expansion [cf. Eq. (1) of Ref.[3. The expansion coefficients An are allowed to be complex. 
Accepting the DG challenge, we consider now the general response instead of the linear limit of I. We do no more 
single out the ground state |5'o), unlike e. g., Eq. (10) of I, and consider the general matrix elements of the relevant 
operators [external bias Hamiltonian, W , and current operator at site q, jq = ietq{a^^aq^i — a^^iaq)] 

Wn„, = (*„ IM^I ^rn) = Km = (*n {a^a^l ^-m) = X;^n. (1) 

(^Ji \ ijq ^ jqo) \ ^m) — ^nm ^ -^nm ~ ^J^mm (2) 



where go is an arbitrary fixed site. (Unless otherwise specified, we use throughout the notations and definitions of I 
and ignore the spin for simplicity.) Notice that all the above calligraphic symbols denote real quantities. To determine 
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the steady-state current, the DG variational ansatz prescribes a constrained energy minimization 
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where w, A^, and Xq ^-re real Lagrange multipliers. The constraints are (see I): 
The minimization with respect to An and ^* of Eq. ([31) yields 
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respectively. Notice the reversed sign of the last term in the square parentheses of Eqs. (O and 
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It is due the 



antisymmetry expressed in Eq. 
operator are purely imaginary, i. 
corresponds to a current + J, j^*' 
(and the Lagrange multipliers), 
computed 



([5|). which is related to the fact that the matrix elements the Hermitian current 
e., I^,^ are real. Rephrasing, this reflects the fact that if the state j^") = J^n ^nl^n) 
) = J2n ^nl^n) corrcspouds to a current — J. To determine the expansion coefficients 
one must solve Eqs. (j4|) — ([8]), and then the position-independent current can be 
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As visible in Eqs. ([7]) and (jH)), An and A* obey different equations. So, the A„'s are allowed to be complex, hence 
1 4*) is also allowed to be complex. Constraining populations, Eq. Q, allows a broken time- reversal symmetry. To 
see whether the ^„'s are indeed complex (thence whether a current J, Eq. (jH]), can indeed flow), one has to solve 
the set of Eqs. dH) — ([SJ. Beyond linear response, Eqs. (0]) — ^ are coupled nonlinear equations, and the solution 
is not necessarily unique. In I, we do solve this problem for linear response, show that the solution is unique, 
and the unphysical result J = is the outcome of these calculations. Throughout our critical analysis of the DG 
approach,i2ili^ to be on the safe side, we confined ourselves to the linear response limit, wherein the solution is unique 
and its lamentable failure can be unambiguously stated. It is possible that the DG current vanishes beyond the linear 
response limit;^ but we cannot safely state this. But if this were the case, it would have nothing to do with the 
DG claim: as seen above, constraining populations docs not prevent a broken time-reversal symmetry. Whether the 
solution of the nonlinear problem corresponds to J = or not did and does not represent our concern once the linear 
response limit is incorrect. This would be yet another unphysical DG prediction, as we already stated explicitly: see 
the last but one paragraph of Sec. VI in Rcf. 0. 

Let us now see what is wrong with the DG argumentation yielding the "conclusion" that constraining populations 
is unphysical. There are two reasons why their conclusion is incorrect: on one hand, the DG argumentation eliminates 
an important aspect and, on the other hand, it introduces an unjustified assumption. DG seem to have realized that 
current conservation is a serious stumbling block for their argumentation; noteworthy, the steps 1 — 6 listed in the 
beginning of the Comment omit the current conservation. They first discuss the minimization ignoring the current 
conservation. This amounts to exclude the last term of the square parentheses of Eqs. ([3]), ([7]), and ([5]) above. This 
almost sufRccs to argue that Eqs. Q — ((8)) automatically imply real coefScicnts An- "Almost" but not really, because 
the unknown An can be complex (Re^„ x Im^„ ^ 0) despite the fact that all the other quantities entering these 
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equations are real (remember x'^ + x + 1 = 0). Then, to more "convincingly" argue that our constraints do not 
break the time-reversal symmetry, DG supplement their argumentation with an unjustified supposition. Namely, they 
simply postulate that "the variational approach assumes a unique solution 4' ..." The DG variational approach is a 
well-defined mathematical problem. Needless to say, whether a solution exists and is furthermore unique cannot be 
"assumed" in any mathematical problem. This should be demonstrated, an imperious requirement particularly in the 
case of such an approach, whose shortcomings are so amply documented.— i^"^ In the last part of their argumentation, 
DG claim that their deduction that ^ is real remains true if they impose current conservation. As discussed above, 
this unjustified assertion is contradicted by our Eqs. ([7]) and ([S]). As a matter of fact, it is just the opposite sign 
of the last term in the square parentheses of Eqs. ([7]) and ((8)) that is related to the broken time-reversal symmetry 
expressed by the transformation 'I' — >■ 'I'*, J — > —J {I -> — / in DG's notation) noted above. DG also mention 
this transformation but fail (or avoid) to correctly include the current conservation as a constraint in the energy 
minimzation that has to be done. 

To summarize, constraining populations is very possible and does not in the least exclude a broken time-reversal 
symmetry. Were our choice to constrain electronic populations inadequate, which are used in the most well-established 
approaches for the uncorrelated and correlated transport [like those based on Boltzmann's equations, nonequilibrium 
Green's functions (NEGF), master equations, etc], all these approaches would fail, but this is not the case. They all 
employ the Fermi distributions (FDs) = (^'In^l^') to express open boundary conditions for incoming electrons. 

Before ending this section, we still make the following remark. Most important for the invalidity of the variational 
DG approach demonstrated in I is that it leads (even with the constraints most justifiable physically, namely n^) to 
an unphysical result. We paid and pay little attention to the fact that the deduced current vanishes. By chance, the 
current could have been nonvanishing, as could fortuitously (but must not) happen when using the DG variational 
ansatz with Wigner constraints (e. g.. Fig. 3 of Ref. |^. There are certainly many possibilities to obtain nonvanishing 
currents by constraining ad hoc complex Hermitian operators, since all these generally yield J 7^ 0: see Eqs. (18)- 
(20) of I for 3^K,Tt 7^ 0; one needs not be too "careful" for this, contrary to the impression which DG attempt to 
convey. Essential to deduce results which are physically relevant is to be careful in selecting from the very broad 
class C mentioned in the Comment those Hermitian operators which are associated to observables able to express the 
physical reality at the boundaries. And, according to the present community's wisdom, a choice better justifiable 
physically than the electronic momentum distributions (FDs) is not known; DG themselves have to admit that they 
represent "an interesting alternative to the" WF (cf. last paragraph of the Comment). 

III. REFUTING THE CLAIM THAT WE USED REAL ORBITALS 

Attempting to give a physical basis to their criticism, DG switch to the single-particle description and claim that 
we cannot constrain only incoming electrons, as appropriate for transport, because our single-particle wave functions 
(orbitals) are real, and outgoing electrons become concomitantly constrained. 

In I (as well as in Sec. |lT]) we need not orbitals at all, nor constrain orbitals, so it is not true what DG write that 
we "have decided to constrain the occupation of real states of the left electrodes of approximate form sin(fc7ra;/L)" . 
These (approximate) wave functions pertain to an isolated electrode (i. e., uncoupled to a device) with open boundary 
conditions (oBCs). Here, open is meant not in the sense of transport theories (of a system exchanging particles 
and energy with environment), but in that used, e. g., in numerical exact diagonalization or DMRG (density matrix 
renormalization group) studies, of a system with open ends whose wave function vanishes at the boundaries [in DG's 
notation, a^^{x = 0) = a^^{x = L) = 0]. Our boundary conditions (BCs) for the isolated electrodes [cf. Eq. (6) of 
I] are general (see Sec. ITV)) . 

Nevertheless, to understand why the DG claim is incorrect let us suppose in this section that we would have imposed 
oBCs. The DG's analysis related to the functions a^^ (x) explains nothing but the trivial fact that a current cannot 
flow in an isolated electrode. How could a current flow in an electrode uncoupled to device? What we constrained 
is clearly visible in Eq. (23) of I: that the populations of the incoming electrons in two many-body states (|^'o) for 
V ^ and [vP) for V 0) pertaining to the total system (electrodes coupled to the device) are equal. We did not 
constrain noninteracting electrons confined within an isolated electrode described by the wave function denoted by 
DG by (x) « sin(27rK/L). This would be meaningless physically. That our constraints refer to incoming electrons 
was already indicated in the step (ii) of Sec. II of I. Our Erratum,— which DG attempt to misinterpret to gain more 
credibility, merely emphasizes this fact; it neither corrects an error nor retracts a statement. 

In fact, as already noted above and discussed in detail below, our results apply to more general cases than that 
analyzed in the Comment, but let us still remain in the latter framework. DG calculations would also be possible by 
using single-electron wave functions (j)i^{x;V) to build all Slater determinants needed to exhaus the total electrode- 
device Hilbert space and express the general many-body state ^' required for the prescribed minimization. (Notice 
that our demonstration of I is exact and does not rely upon any approximate CI expansion, contrary to Ref. |3|) . The 
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physical content of the orbitals (t>n{x, V = 0) with (t>^{x = qL:V = 0) = Q [these are the DCs {x)] is that incoming 
(in) and outgoing [out) electrons move symmetrically in the isolated left electrode. If we formulated in I the DG 
variational ansatz in terms of orbitals, the transport approach would have allowed a broken m(coming)-oMt(going)- 
symmetry when the electrodes are coupled to a device o„i(a; = qL',V) 7^ 0] for a nonvanishing bias V ^ Q. K 
current J can flow only through a device coupled to electrodes. These (ji^ix', T^)'s would have been used throughout 
to express all averages needed in I [including the rhs of Eq. (23), that is, at = for electrodes coupled to a device] 
and determined self-consistently from the DG variational ansatz, just like our An& of I or above in Sec. |TT1 Whether 
this symmetry, allowed to be broken, is indeed broken and a current indeed flows (J 7^ 0) should be determined 
by performing calculations within the DG variational ansatz, and this is exactly what we did within the many-body 
formalism of I. // the DG variational ansatz were valid, the minimization would yield wave functions (x; V ^ 0) 
different from (p^aut {^'^ ^ 7^ 0) ^^d J ^ 0, in the same way in which it would have led to ImA„ 7^ in I and in Sec. [IT] 
above. Of course, this cumbersome approach based on the first quantization is equivalent to the second quantization 
formalism employed in I, and what is important is the unphysical outcome (J = 0). It clearly demonstrates that the 
DG transport ansatz is invalid. 

IV. REBUTTING THE CLAIM THAT COMPLEX ORBITALS RESCUE THE DG VARIATIONAL 

ANSATZ 

DG chose the oBCs as a stumbling block for the demonstration of I, but, because that case turned out to be more 
subtle, ended becoming trapped themselves in this pitfall. By contrast, the case of pBCs indicated by them as a 
way-out solution is really much easier to analyze and allows to easily understand the correctness of the result J — Q 



As already noted, our derivation of the unphysical result J = for the linear response limit in I is very general. To 
give a flavor on how general are the situations where the DG variational ansatz fails, we referred to a rather broad 
class of uncorrelated and correlated models, which are mostly used in nanotransport, although our demonstration 
applies well beyond that class. In I we gave a formal general analytical demonstration and needed not bother to do 
expensive numerical calculations for certain device models nor become involved by giving a particular type of the BCs 
for the isolated electrodes or explicitly indicating which labels k refer to single-particle states of incoming electrons 
and which to outgoing electrons, although they can be obviously specified (see below in this section and Sec. IVip . 

By inspecting the second and third lines of Eq. (6) of I, one can immediately realize that our isolated electrodes can 
be described by general boundary conditions. Just for not to bother the reader with an unnecessary lengthy discussion 
of the BCs, we avoided, e. g., to specify the site indices at the ends of the electrodes opposite to the device, and 
wrote / < qL and r > g/j in the second and third lines of Eq. (6) of I. The message of our Refs. is unambiguous 
and leaves no hope: the DG approach is completely incorrect and cannot be rescued. Having demonstrated this, we 
doubted and doubt on the usefulness or necessity of a too detailed analysis of an approach so incorrect. 

The results of I apply for the aforementioned oBCs (Sec. IIIip . but also for boundaries that are, e. g., periodic 
(pBCs), anti-periodic (Mocbius) or general twisted (see, e. g., Ref. |^ and citations therein). To this, one should 
keep, e. g., Mn sites in the right electrode and impose, e. g., that the annihilation operators satisfy Oq^^^Ma = o-qn, 
aqj^+Mn = ~o-qm or o-qa+Ma = oiqu exp(if?) {6 is a real constant phase), respectively. The left and right electrodes can 
even be described by different types of boundaries. 

In their Comment, DG concede that our idea of constraining populations "is an interesting alternative to the Wigner 
function" , but claim that it is essential to choose complex single-particle wave functions associated to pBCs. Accepting 
also this DCs challenge now, let us explicitly work out in detail just this case where the pBCs arc applied to the isolated 
electrodes and where the orbitals are complex, which DG suggest as a way out to rescue their variational ansatz. From 
now on, we shall always assume pBCs unless otherwise specified. With pBCs, and assuming homogeneous hopping 
integrals Tl,r in electrodes for simplicity, one can write down explicit analytical formulas instead of the general 
unspecified ones of I, because the transformation matrix T which diagonalizes H^ n entering Eq. (22) of I is nothing 
but the well-known Fourier transformation 



of I. 



O-qL+l = 



,2TTikl/ML 



k 




(10) 



p 



Throughout, we use I = -Ml + 1, . . . , 1, 0; r = 0, l,...,Mfl - 1; -Ml/2 < k < Ml/2; -Mr/2 < p < Mr/2, 
and assume even Ml.r for specificity. With Eq. pUj) . the Hamiltonians Hl, Hr, and Hd^c of Eq. (6) of I become 
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(tl = tq^ , and tr = tq„ ) 

Hl = ^[fiL - 2TLCOs{2Trk/ML)]alak, 

k 

Hr = ^ [ma - 27r cos{2itp/Mr)] alup, 
p 

fe 
p 

Further operators affected by the transformation pop entering the relevant equations of I are 

fe <3 p 

fe 

jR = i|rLA/^^/^^ (a],a,^ - aj^ap) , (14) 
p 

representing the Hamiltonian of the external bias W and the current operators jL,R at the contacts, respectively. The 
sites within the device g = + 1, . . . , g/? — 1 are not affected by Eq. ([TU|) . and the corresponding operators (device's 
Hamiltonian Ho and the current within the device jq) can be found in I. 

By inspecting Eq. ([TT]) . one can immediately see that, in spite of the fact that the single-particle wave functions 
i'kixi) ^ exp(27rifcZ/Afi) and (j)p(xr) ~ exp{2T:ipr / Mr) of Eq. (|10p are complex and degenerate {+k and —k correspond 
to the same energy), all the parameters of the Hamiltonian H = Hl + Hr + Ho,e + Hd are real. Consequently, 
all its many-body eigenstates jvPn) are (can be chosen) real. All the parameters entering W are real, so the matrix 
elements (*„|VF|*o) = W„ [Eq. (13) of Ref. are again real. As in the general case of I, the matrix elements of 
the current operator are purely imaginary (i. e, real JTq^n)- AH the matrix elements of the particle number operators 
ajla^ (k = p, fc) remain real (4'„|a]^aK|^o) = -^k (3^k = 0). Eqs. (24) and (25) of I remain unaltered. Whether 
imposing current conservation (as we did) or not (as DG incorrectly claim that one could do^, as if their method were 
so good to automatically include current conservation^), the completely unphysical result (J = 0) follows as the 
ineluctable conclusion of applying the DG variational ansatz. For these pBCs, the labels of the single-particle states 
of incoming electrons can be explicitly given (see also Sec. IVip : — >■ at^^ with < km < Ml/2 and — >■ ap.„ 
with —Mr/2 < Pin < 0. One can now explicitly see that only these incoming electrons can and are to be constrained 
in Eqs. (12) and (23) of I or in the present Eq. They do differ from the outgoing electrons, whose labels are 
—Ml/2 < kout < and < Pout < Mr/2, and one can convince oneself explicitly that outgoing electrons (can) 
remain unconstrained. So, these constraints correspond to Fig. 7c of Ref. and not to Fig. 7a and 7b, contrary 
to what the Comment claims. By simple algebraic manipulations of Eqs. (23)-(25) and (15) of I one can easily 
deduce that the distributions of the outgoing and incoming electrons are equal, (^'|aj^afc|4') = (vtlal^^a-fcl^) and 

(4'|aLpQ:_p|\I') = (\E'|aJap|4'). (Of course, this is not at all surprising in view of the unphysical result J = 0.) One 
can now explicitly see that the equality of these distributions is the outcome of the calculations of the defective DG 
variational ansatz, and is by no means (not even implicitly) assumed from the very beginning through an inappropriate 
choice. The constraints used by us within the calculations based on the DG variational ansatz allowed a broken 
symmetry between incoming and outgoing electrons. Whether this symmetry is broken or not remained an open 
result, which emerged from the DG transport calculations. The result is that this symmetry is not broken, J — 0, 
demonstrating the incorrectness of the DG variational ansatz, and this also holds true for the pBCs, contrary to what 
the Comment argues. From the above analysis of the pBCs it is also clear that the electrodes' size Ml,r, can be 
arbitrary large (which is impossible within the ab initio DG calculations^). Therefore one can also understand that 
the unphysical prediction J = is not limited to pBCs but holds for any other BCs as well, since otherwise, e. g., the 
entire philosophy of solid-state physics to apply pBCs would break down. 



V. WIGNER FUNCTION CONSTRAINTS VERSUS PARTICLE DISTRIBUTION CONSTRAINTS 



DG mention (as we also did) that, by constraining the WF, a nonzero current is possible. In Refs. [^HS we 
discussed that, luckily, mathematically this may be possible. Why did DG constrain the WF? Only because in Ref. H 
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they claimed that the FDs cannot be used for correlated transport. This claim does not at least apply for uncorrelated 
systems. At least there, FD-constraints are possible, and our critique of the DG variational ansatz from I obviously 
applies. In fact, in I we explained that open boundary conditions can also be formulated by means of FDs even for 
correlated electronic devices, because they should be imposed in electrodes, wherein electrons are uncorrelated. It is 
the FD which has a precise physical meaning, and not the WF, which has a physical content only when it is a good 
approximation for the FD. According to the Comment's original philosophy, a current flow is possible only to the 
extent to which the WF does differ from the FD. 

The inappropriateness of the Wigner constraints could not be immediately recognized only because, luckily, the 
matrix elements of the Fano operator are generally complex. They are generally complex no matter whether the 
single-electron functions are real or complex (see Refs. 011); this is not the result of any "careful" choice of certain 
complex functions, as incorrectly claimed by DG. Let us show that this is also the case when the Comment's continuous 
space description is used instead of the discrete one of Refs. [^Q- The Fano operator reads 

F{x,p) = ^jdr e-'P^'i^^x - r/2)i;{x + r/2), (15) 

where ■!/'^(x) and '4>{x) the electron field creation and destruction operators. Its general matrix element for two arbitrary 
A^-body states \^) and |$) corresponding to the multielectronic wave functions \l/(a;i, . . . , xat) = {xi...xn\^) = 
{Q\ij^{xi) . . . '!/>(xAr)|^')/-\/iVT (|0) is the vacuum) and $(a;i, . . . , xn) = {xi . . . xn\^) can be easily expressed as (Xjv-i = 

{cCi, . . .,XN~l}) 

X *(xAr_i,.T + 0dXw-idr. (16) 

This matrix element is generally complex irrespective of whether the wave functions 5'(XAr) and $(XAr) entering 
Eq. (|16p arc real or complex. Eq. (jl6p is general and holds whatever the employed single-particle wave functions 
(which need not be specified there). The complex exponential entering Eq. (|16|) belongs to the definition of the Fano 
operator, Eq. ([T5|). and has nothing to do with the employed orbitals. Whatever the latter, it can be artificially split 
as exp(— ipr) = exp[ip(a; — r/2)] x exp[— ip(a; -I- f /2)]. This trivial splitting and the notation a^'^(a;) for these factors 
is obviously done by DG in their Eq. (7) merely for conveying the false impression that complex exponentials would 
represent a key point, which they would have "carefully" exploited in their work.— 

Let us express the population constraints for the incoming electrons [p < 0) of the right electrode (TZ), Eq. (23) of 
I, making use of the Fano operator— and Eqs. ([TU]) 

^{^\F{x,p)\^) = (^\alap\^) = {^o\alap\^o) = ^{^o\Fix,p)\^o) ■ (17) 

Instead of constraining the above sums (integrals) over x of the Fano operators (populations), DG's discretionary 
constraint of a single term in either electrode (namely, x = qL,Ft.)^ involves a quantity which does not possess a 
physical meaning. Emerging from such an ad hoc mathematical constraint, it is not at all surprising that the currents 
predicted by the original DG approach;^ whether they vanish or not, arc completely unphysical, as demonstrated in 
Refs. Hi. 

To conclude, mathematically the WF-constraints used in conjunction with the DG variational ansatz (can) yield 
a nonvanishing J (without any physical relevance, cf. Refs. 0111) even when using real wave functions, in spite of the 
claimed "warning bells" that "real wavefunctions carry no current" »i 



VI. CONSTRAINTS FOR n-TYPE AND p-TYPE CONDUCTION 

In Sec. IIVI to specify the labels of incoming and outgoing electrons, we have made the intuitive (or, better, naive) 
assumption that the single-particle states with positive (negative) wave numbers k and p correspond to right- (left-) 
motion. However, they are related via Eqs. (jlOp to gwasz-momenta and not necessarily to physical momenta. Especially 
for later purposes, it is important to demonstrate that their sign is indeed related to the direction of the motion in 
the real world. 

Let us consider the isolated left electrode. Its ground state (not to be confounded with that of the coupled 
electrode-device system, |4'o), cf. Sec.lIIip is the Fermi sea \F) = [Y\iu<^k^ |0), where the Fermi wave vector kp is 
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determined by the number of electrons. Using Eq. (|10p one can straightfordwardly demonstrate that at any position 
(iL + 1 within the (left) electrode, the average of the electron number current jq^+i = *(agj^_|_/_|.ia<ji+; — h.c.) vanishes, 
(^\^qL+i \^) ~ ^- [Notice the opposite signs of j and the electric current j of Sec. Illfor electrons.] Let us also consider 

the states [Ml/2 > \K\ > kp, \K'\ < kp < Ml/2) 

\<f^^)^aUF),\<i>'k,)^aK,\F). (18) 

They represent states with one extra electron (el) and hole (h) in the Fermi sea, respectively. Straightforward 
calculations using Eqs. ((T8| and ((TO)) yield 

wiJ„+,i*i;> = +2^»i"^. (19) 



The sign of the j-average does express the real direction of electron quantum-mechanical motion. Therefore, Eqs. (fT9|) 
and (PU]) demonstrate that the sign of the wave vectors belonging to the Brillouin zones of Sec. IIVI (symmetric around 
zero) specifies the direction of electron motion, and that electrons and holes with a given wave vector move in opposite 
directions. The latter result can also be seen by performing the general particle- hole transformation, e — > — e (charge 
conjugation) and {ip^x),ip{x)} {ipl{x) = 'il!{x),ijjh{x) = ip^ix)} {{al,a^} {a'^^ = ax,a'^ = aj,}). Using Eq. ([15]), 
one easily gets 

jt ^ -jx'Jx ^ (21) 
Fh{x,-P) = ^ j dr e'P''ii{x~r/2)Mx + r/2) = -Fix, P)+ const. (22) 

In view of the aforementioned, one can conclude that incoming and outgoing electrons correspond to the wave vectors 

' I'To^ ^ ';o(for electrons) 

-Ml/2 < kout < 0; < Pout < Mr/2} ' (24) 



while for incoming and outgoing lioles 



-Ml/2 < fc,„ < 0; < < Mr/2, , , 
< kout < Ml/2; -Mr/2 < p^ut < 0> ' 



(25) 
(26) 



The fact that the above electron and hole descriptions are equivalent is trivial in general, but not in the context of 
transport approaches, wherein incoming cfiarge carriers arc to be constrained.— If the charge carriers arc electrons 
(n-type conduction), the constraints should be imposed to incoming electrons, Eq. (|23p . In this case, the intuitive 
assumption of Sec. IIVI is justified. However, if the charge carriers are holes (p-type conduction), one should constrain 
the incoming holes, Eq. ([25]); that is, the labels k in the above Eq. ([4]) and in Eqs. (12) and (23) of I are those given 
by Eq. ([25]) and not by Eq. ([23]). 

The analysis of this section and of Sees . IIIII and IIVI makes it now clear why we preferred to consider the general case 
in I and not to enter in unnecessary involved details: they are absolutely not necessary to understand the unphysical 
prediction J = of the DG variational approach, and hence its lamentable failure. But because the incorrect DG 
claims in the Comment brought us to enter such details, we can show another shortcoming of the original DG approach 
related to them, which we did not present so far. 

In their workj^ DG did not examine at all whether the molecule they considered, BDT (benzenedithiolate), exhibits 
an n- or a p-type conduction. Uncritically, they merely constrained f{qL,P > 0) and f{qR,P < 0). Even if their 
variational ansatz were correct, and even if these WFs were true distribution functions, these constraints would be 
appropriate only if the charge carriers were electrons [conduction mediated by LUMO (lowest unoccupied molecular 
orbital)]. In reality, in BDT the majority charge carriers are holes (p-type conduction), as clearly demonstrated by 
the recent, accurate experiment of Ref. [lOj- By inspecting now Eqs. ([22]) (noting the reversed sign of P in the Ihs 
and rhs), (|23p . and ([26]), one is amazed to see that what DG constrained in Ref. |3| are in fact the outgoing majority 
carriers, and not the incoming ones. It is certainly too simplistic to describe the conduction through BDT merely 
as a process mediated by HOMO (highest occupied molecular orbital; p-type conduction) instead of accounting for 
several/numerous ionization and electroaffinity levels, but the fact that the constraints of majority carriers (holes) 
are unphysical in Ref. is a clear demonstration that uncritically using Wigner boundaries is completely unjustified. 
To conclude, even if all the other DG ingredients were correct (what is obviously not the case^), this very reason 
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irrefutably demonstrates that the results of Refs. I3I IIIII12I have absolutely no physical meaning. What would be the 
appropriate constraints in the case of ambipolar conduction, where both electrons and holes contribute to the current, 
is an issue,— which we do not discuss here. 

In our first work^ that challenged the DG approach with Wigner constraints;^ we considered uncorrelated and 
correlated quantum dots modeled by a single level whose energy offset from the electrodes' Fermi level sp = y^L = fJ-R 
is Eg. The results of the DG calculations presented there (e. g., in Figs. 2-5 and 7) are for Eg > 0, that is, the dot's level 
plays the role of a LUMO (n-type conduction). The charge carriers are electrons, and our constraints [corresponding 
to the above Eq. ([25)1 ] refer to incoming electrons. 

Both the uncorrelated and the correlated models of Refjj are described by Hamiltonians H{Eg) possessing a particle- 
hole (or charge conjugation) symmetry (see, e. g., Ref. [l3l and citations therein) around Eg = 0: Hh{Eg) = eA 
That is, the zero-bias conductance g[-\-£g) = g{—£g) (as well as other relevant properties not considered in Ref. 13, 
e. g., the whole current-voltage characteristics) should be identical irrespective whether the level is located above 
i+Eg) or below {—Eg) the electrodes' Fermi level Ep- Noteworthy, the charge carriers are electrons for positive Eg 
and holes for negative Eg. As a test for numerical calculations, we checked that DG calculations for the LUMO case 
{Eg > 0) constraining the incoming electrons [Eq. (|23|) ] and for the HOMO case (e^ < 0) constraining the incoming 
holes [Eq. (j25|) ] yield the same, albeit completely unphysical linear conductance. [The electric current operator has 
the same sign both in the electron and the hole representation, cf. Eq. (|2ip .] As clearly demonstrated,— the DG- 
conductance computed in this way is completely unphysical, but . . .it is still positive, gDG{^g) ^ both for positive 
and negative Eg. That is, this (modified) DG approach can still "predict" that electrons flow from the lower potential 
to the higher potential, and holes flow from the higher potential to the lower potential. 

If we drew the curves of Figs. 3, 5, and 7 of Ref. also for Eg < 0^ by blindly computing the DG conductance using 
exactly the DG prescribed constraints^ [i. e., Eq. ([13])] we could have shown a funny "prediction" of the DG approach^, 
namely, that the linear conductance can be negative, gDG{£g < 0) < 0! That is, holes should have to flow . . .from 
the lower potential to the higher potential. This results from the fact that the blind constraints of f{qL,P > 0) and 
f{qR,P < 0) erroneously constrain in fact the outgoing carriers; this situation corresponds to Fig. 7d, and not to 
Fig. 7c of Ref. i. Indeed, these DG Wigner constraints break the time-reversal symmetry and yield a nonvanishing 
current, but . . . what is the physical relevance? As a matter of fact, it is just such an unphysical imbalance, which 
is shown in Fig. 1 (bottom) of Ref. (the counterpart of Fig. 7d of Ref. |^ and not of Fig. 7c, as incorrectly claimed 
in the Gomment), that breaks the time- reversal symmetry in Ref. In Ref. i, we did not show this conductance 
g{Eg < 0) < because the demonstration of the severe failure of the DG approach with Wigner constraints was 
sufficiently convincing even without mentioning this "prediction" , and we preferred to couch the discussion in terms 
as sober as possible. However, we have noted it above, since, in spite of the clear evidence of Refs. 011, DG still 
continue to uncritically refer to their work^ in the Comment. 

VII. FURTHER ERRORS AND INACCURACIES IN THE COMMENT 

In the Gomment, DG claim that we "doubt the validity of using the Wigner function constraints to apply open 
system boundary conditions" . As we repeatedly emphasized^i^i^ we did not challenge in any of our works published so 
fari^jli^ the imposition of Wigner constraints in general; we irrefutably demonstrated that applying Wigner constraints 
in the specific context of the DG variational approach yields completely incorrect results. 

Another inaccurate assertion of DG is that "Baldea and Koppel accept that our [i. e., DCs] Wigner constraints 
. . . have . . . non-zero linear response current." In reality, we demonstrated (i) that even if, by chance, this is possible^ 
the DG current is completely unphysical, and (ii) that the DG conductance vanishes just in on-resonance cases, where 
the exact conductance attains its maximum (unitary limit) 

To refute the false impression, which DG attempt to convey by using the superscript DG, we emphasize that they 
did never impose pBCs for the Awia-clusters, which mimic their electrodes of Refs. MTllfl^ . For this, it suffices to 
inspect Eqs. (I) and (2) of Ref. i. Still, let us emphasize for completeness that, even if they imposed such Wigner 
pBGs, the results of the DG approach would have been incorrect: in Sec. V of Ref. we demonstrated that the DG 
conductance computed with Wigner constraints vanishes {goG = 0) also with pBCs for electrodes, and this occurs just 
in a typical, well-known physical situation corresponding to the unitary limit, wherein the true conductance reaches 
the maximum value {gexact = e'^/h). 

Our Ref. H is the only publication wherein the DG approach, as originally proposed (i. e., with WF-constraints) , was 
worked out using plane wave orbitals exp(iKx) in electrodes with pBCs. DG used nowhere the plane waves cxp{iKx) 
mentioned in the Comment or other orbitals pertaining to electrodes with pBCs, which they denote by a^^{x) in 
the Comment. Or, more precisely, they should have denoted them so [i. e., aj^'^{x) = {x\a^'^) = exp{iKx)], because 
in fact their expressions, like [a^*^) = exp{iKx) in the last but one paragraph of the Comment, do not comply with 
textbooks' quantum mechanics. Concerning the operators denoted by Q^'-^ , they are useless: these operators are not 
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defined at all. 

Most significantly, the Comment contains equations and assertions, which defy more than eight decades of using 
the second quantization formalism. E. g., DG incorrectly state that the action of the operator Qf^^ = '^l^^'^k^ 
on a one-electron wave function ipix) is "straight- forward" . This operator cannot act at all on ip{x), not even on a 
many-electron wave function X2, . . . , xn) in the coordinate space, as expressed in their incorrect Eqs. (1) and (2), 

respectively. It acts on many-body states belonging to the abstract Fock vector space. Further, it is incorrect what 
they write that the operators aj.^^ and aj^^ create and destroy the ei^erastates of ■ Likewise, these operators act 
on any many-body Fock states {'^general) , wherein they create and destroy one electron in the single-particle state k. 
We thank DG for using the superscript BK in their notation al^^ and a^^ , but it is superfluous. The creation and 
annihilation operators a' and a used in I [the same as those entering the present Eqs. (jlOp ] are already unequivocally 
defined in textbooks, and it is in this sense that we used them throughout. 

In view of the aforementioned, the need to translate in the Comment obvious results expressed in the second 
quantization in I into the first quantization language [e. g., the attempt to "demonstrate" that the matrix element 
{'^n\Q'!^^ l^^o) of Eq. (4) in the Comment is real] is understandable. But then it is not surprising that understanding 
what are the actual constraints of I or the generality of the demonstration of I is not obvious. Remarkably, although 
DG reproduce our expression of in the Comment, they completely overlook the content of this expression: namely, 
the fact the BCs are general, and so are the pertaining orbitals, which are not necessarily real. Not coincidental is 
also the fact that DG indicate as a way-out solution just a particular case of I, wherein the incorrectness of the DG 
variational ansatz can be immediately understood (cf. Sec. IIV[) . 

Above, we preferred to respond to and rebut in detail all the issued raised by the Comment. In fact, this represents 
the main objective of our Reply. Still, we note that this analysis of all the concrete aspects is actually not necessary. 
In Sec. VII of our recent work^ we pointed out more serious reasons why the DG variational approach fails. This 
criticism also applies to the case discussed above; it comprises fundamental aspects not confined to a certain type of 
BCs, let they be in terms of WFs, FDs or others. Our criticism was presented in Ref. [^in sufficient detail and will not 
be even summarized in this Reply. A further fundamental reason why the DG approach is incorrect will be presented 
in Sec. El 



VIII. THE WIGNER FUNCTION IS UNSUITABLE TO SPECIFY THE DIRECTION OF MOTION 

The incorrect claims of DG made us aware of a limitation of the usefulness of the WF for transport, which we could 
not find in the literature. 

The WF (/) is employed in many physical studies, including transport's, in spite of its physical limitations. The 
limitation known from textbook a^^i^^ traces back to the Heisenberg's uncertainty principle. The WF can be negative 
and should be not interpreted as a probability distribution, but rather as "one step in the calculation . . . never the 
last step, since" is not measurable but "is used to calculate other quantities that can be measured . . . the particle 
density and current" and "no problems are encountered as long as one avoids interpreting / as a probability density" 
(quotations from ch. 3.8, p. 203 of Ref . [l^ . 

As noted above, in transport it is helpful to distinguish between incoming and outgoing electrons. To this aim, it 
is necessary to use a physical property enabling to indubitably assess that electrons are, say, left- or right-moving. 
The averages of the particle current operator used above or the physical momentum do represent such properties. 
To see whether the momentum "variable" P of the WF f{x, P) justifies to speak of left- or right-moving electrons 
depending on the sign of P, let us consider N noninteracting electrons confined within a one-dimensional square well 
of width L and infinite height. (This could be the isolated electrode considered by DG.) Electrons occupy energy 
levels h^K^ /{2m), whose single-electron wave functions 4>k{x) = (2/ L)^^^ sui{kx) = i(2L)"^/^[exp(— ikx) — exp(zKx)] 
(k — > Kn = Tm/L, n ~ 1, 2, 3, . . .) are just those that DG^ incorrectly claim we would have used in Ref. d. In the 
ground state, the lowest N levels are occupied up to the Fermi "momentum" pp = hkp = nhN/L. Computing the 
Wigner function of this system is straightforward: f{x, P) = X]k„<kf /2-±r/2<o ^'''^■'^P(^*^'')'^'^n {x + r/2) 

(see, e. g., Ref. [ll ch. 3.7, pp. 202-203). 

One might think that one could use the WF as if it were a distribution function in cases where its shape resembles 
a Fermi distribution. Let us inspect the curves for /(x, P) computed as indicated above and presented in Fig. [TJ In 
fact, at smaller sizes (close to the linear size of the DG's Auia-clusters^) the WF does not bear much resemblance 
to a Fermi function (the lower curves of Fig. [T]), At larger sizes (much larger than those one could hope to tackle 
within ab initio calculations to correlated molecules, for which the DG approach^ was conceived) the curves (the upper 
part of Fig. [T]) become more similar to a step function, and one may think that this is encouraging. In reality, the 
contrary is true: as visible in Fig.[TJ mathematically one can calculate the WF for positive and negative "momentum" 
variables P separately. However, this mathematical separation does not reflect a physical reality: for any single-particle 
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FIG. 1: (Color online) Wigner function for two sizes L computed at two points x indicated in the legend. Gold's Fermi wave 
vector pp/h— 12nm~^ is used. 



eigenstate k the electron momentum vanishes, = —ih J dx(/)%{x) (d/dx) 4>k,{x) = 0; left- and right-traveling waves 
are entangled with equal weight, and one cannot speak of single-particle eigenstates representing left- or right-moving 
electrons only because Wigner functions with positive or negative P-argumcnts can be computed. This represents a 
further limitation of the usefulness of the WF, not related to the Heisenberg's principle, which is particularly relevant 
for transport. The current has a direction, and if one wants to unambiguously specify this direction, the WF f{x, P) is 
inappropriate; a WF with negative (positive) "momentum", P < (P > 0), does not imply that left- and right- moving 
particles exist in the real physical world. 

So, using /{ql, P > 0) and f{qR, P < 0) as if they were true momentum distributions of incoming electrons, as DG 
did,— is not justified in quantum mechanics. The above example demonstrates that, indeed, the textbook's warning 
mentioned in the beginning of this section is pertinent. 



IX. WHY ANY APPROACH TO TRANSPORT MERELY BASED ON A FINITE ISOLATED CLUSTER 

NECESSARILY FAILS 

In the course of our extensive critical investigations^ii^"— of the DG approach^ we became aware of a series of 
difficulties, which not only the DG's but also other approaches to nanotransport are faced with, which we want to 
bring to the reader's attention. 

To understand the problem, let us briefly consider the uncorrelated dot model of Ref. linked to semi-infinite 
electrodes {Mi^ n — > oo). Exact transport calculations at arbitrary bias V can be easily carried out within a multitude 
of approaches'^"— Besides the current J and dot occupancy Nq^— the occupancies of the sites in electrodes can 
also be computed 

In Fig. [21 we present steady-state results for the electron number on the dot A^o and that for nanoclusters Nc 
centered on the dot and including several electrodes' sites Ns- As visible there, the changes in both V (source-drain 
voltage) and Eg (gate potential) yield variations in the dot occupancy Nq. (Notice that there is no exchange of 
electrons with the gate.) The total number Nc of electrons in the nanoclusters also varies, it closely follows the change 
of iVo; the variations in the small difference Nc — Nq could be hardly seen within the drawing accuracy of Fig. [2J and 
therefore are not shown there. The fact that Nc — Nq is practically constant implies that the sites in electrodes remain 
practically unaffected by changes in V and e^; even the electrodes' sites in the very close dot's vicinity are very little 
affected. The dependence Nc = Nc{V,eg) is most important in the context of a transport approach. It demonstrates 
that, by varying V and/or Eg, a finite nanocluster exchanges electrons with the infinite electrodes linked to it, which 
act as reservoirs that can supply /withdraw electrons. 

In the light of the aforementioned, it becomes clear that it is strictly impossible to correctly describe the transport 
within a DG-like approach, which considers a small isolated cluster. Within an approach like that initiated by DG^ 
and further scrutinized by us^^i^ one attempts to describe the transport using such a cluster characterized by a 
many-electron wave function 'i!(xi,X2, ■ ■ ■ ,X]\[) to be determined by a constrained minimization of cluster's energy 
for V 0. This wave function merely describes a nanocluster with a given number of electrons N , which cannot be 
varied by changing source-drain (V) or gate (Eg) voltages: 
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FIG. 2: (Color online) Dot charge A'o and cluster excess charge SN^i = Nc — Ns for two gate potentials Eg = 0.5 {Ns = 3, 9, 11) 
and Eg = 0.25 {Ns = 3, 11, 13). For both gate potentials, the employed cluster sizes A'^a given in parentheses increase upwards 
for, say, 1/ ~ 3. Notice that for larger sizes the excess charge rapidly saturates and closely follows the dot charge, and this 
gives a flavor on how fast the convergence to charge neutrality within more distant parts of electrodes is (see the main text). 
Energy unit is Tl,r = 1 and tl,r = 0.2. 

(a) whatever the property chosen to be constrained (WF, electron momentum distribution or else), 

(b) whatever the boundary conditions (oBC, pBCs, Moebius, general twisted or else) which might be chosen for the 
finite "electrodes" included in the cluster (the Auia-clusters of Ref. 0), 

(c) let the many-electron wave function "^{xi, X2, ■ ■ ■ , xn) be determined within exact full CI expansions (i. e, 
using the whole Hilbert space) as done by us^^i^ and not only approximately, by using a reduced number of 
configurations selected by Monte-Carlo sampling, as DG did<^ 

Even the most "carefully" chosen constraints can at most, e. g., acceptably account for the change the number of 
electrons in the device (dot). But this change will inevitably modify the charge of the neighboring sites in electrodes. 
It is easy to imagine how profound will be the impact on the transport, e. g, in correlated or switchable devices. 
It is hard to conceive that the transport could be reasonably described within this framework, even if the cluster 
would be very large (what is obviously hard within the presently available ab initio quantum-chemical calculations for 
correlated systems): because it is by no way obvious/necessary that the local charge density be incorrectly described 
only at the remote ends of the finite "electrodes", which would allow to suppose/hope that the electric field within 
the device will be little affected. A further difficulty is, of course, the fact that the number of electrons Nc{V,eg) of 
the coupled nanocluster is generally noninteger (cf. Fig. [5]), but we do not discuss this issue herei^ 

One may further ask whether a "careful" choice of a certain charge state of the dot (molecule) instead of the 
neutral species could help. No, it does not. For the above model, one can also compute the exact time-dependent dot 
population No{t) by suddenly coupling at i = a dot with population A^o(O) to infinite electrodes.—"— The result for 
y ^ within the wide band limit {4T > £g, F, where Y = 2t'^/T,T = Tl=Tr, and t = tl = tr) ia^ 



As seen in Fig. |3l which visualizes the result expressed by Eq. ((27|) . the asymptotic value of the dot population 
limt_j.oo N{t), which corresponds to the steady state, does not "remember" the initial value A^(0); it is determined by 
the applied voltages. 

The foregoing analysis makes it also clear that not only a variational ansatz like the DCs is incorrect. The above 
list (a) - (c) can be enlarged by stating that, whatever the ansatz (let it be a variational ansatz better than the 
DCs or based on e. g., Liouvillo^ or Schrodinger equations) used to compute the wave function for a fixed N, it will 
fail. The aforementioned flaw is serious and can be unambiguously traced back to the use of a finite (small) isolated 
cluster. Electric transport can only occur in an open system, and allowing electron exchange with environment (infinite 
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FIG. 3: (Color online) Curves for the time-dependent dot charge No{t) showing that the asymptotic values do not keep track 
of the initial occupancy ^(O). The values of the source-drain bias V, energy offset Eg, and hybridisation F are given in the 
legend. 

electrodes) is indispensable. Most commonly, this is done via the embedding self-energies within the Keldysh-NEGF- 
approachfi^ which enable the nanocluster to change the number of electrons it contains. 

X. CONCLUSION 

Most important for the present Reply, we have demonstrated above that the DCs claims of the Comment are 
incorrect and our critique of I is untouched. 

The unphysical current (J = 0) obtained in I is the result of the DG defective ansatz and by no means emerges 
from pretended unsuitable constraints of real orbitals, as DG incorrectly claim. Using complex orbitals does not in 
the least change the unphysical prediction of the DG variational ansatz; furthermore, it confirms more directly its 
incorrectness. In fact, the complex orbitals given by DG represents nothing but a particular case of the general case 
considered in I, for which our demonstration holds. For completeness, we have worked out this case in detail to 
explicitly specify the labels of the incoming charge carriers (not necessarily electrons) , to demonstrate that they can 
be constrained without concomitantly constraining the outgoing carriers, and to show that the distribution of the 
outgoing carrieres represents the outcome of transport calculations. Of course, fully consistent with the unphysical 
prediction J = 0, the distribution of the outgoing carriers deduced within the defective DG ansatz is found equal to 
that of the incoming carriers. 

To conclude, DG concede that our idea of constraining populations "is an interesting alternative to the Wigner 
function", but argue that it is essential for this to choose complex states. We have shown above that the usage of 
the complex orbitals just of the form indicated by them as key point of a way-out solution does not mender the DG 
variational approach. These are not only complex, they are just of the form a^*^(x) = exp^iKx), which DG give 
(and incorrectly attempt to suggest that they used them in Ref. . Consequently, they must now accept that their 
variational ansatz for transport is incorrect. 

The analysis done in conjunction with the issues raised by DG has led us to reveal two aspects of more general 
relevance for the transport theory. First, we have presented an example illustrating a limitation of the Wigner function 
f{x,P) important for transport, namely that the sign of P is not necessarily related to the direction of motion in 
the real world. Second, we have presented an important physical reason why transport approaches, which, like the 
DG's, use information pertaining to a finite isolated cluster are inappropriate. This enlarges the basis of our critique 
recently formulated in Sec. VII of Ref. d, where two other fundamental reasons were exposed. 

To summarize our detailed investigations on the DG approach)^ we can state that this approach lamentably fails 
because virtually all its ingredients are incorrect: 

• The DG approach imposes boundary conditions as if the WF were a true particle distribution, which is not justi- 
fied quantum mechanically. Replacing the WF-constraints by the boundary conditions most justified physically 
(namely, the Fermi distributions) attempted in I does not remedy this approach. 
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• Even if the WF were a true momentum distribution, by uncritically constraining f{qL, P > 0) and f{qR, P < 0) 
it is very possible to erroneously constrain the outgoing majority charge carriers and not incoming ones. (This 
is just the case in Refs.ililli) 

• The DG approach attempts to describe the transport by using a finite cluster within calculations based on a 
variational principle (entropy maximization) , which turned out to be problematic even if, unlike in the DG case, 
the limits of infinite time and infinite volume are taken in the correct order.— 

• The DG approach aims at describing the transport by means of a wave function determined for a (small) isolated 
cluster, whose number of electrons is fixed, while a nanocluster in a real electric circuit does exchange electrons 
with the infinite electrodes to which it is connected. 

The lamentable failure of the DG approach was demonstrated by explicit calculations for the simplest uncorrelatcd 
and correlated, discrete and continuous modelsi^i^ They contradict well-established experimental and theoretical 
results, and it would make little sense to more amply document the incorrectness in many other cases.— By contrast, 
DG could not present even a single example where it is valid. 

Based on ingredients unfounded physically, it is not at all surprising that the currents predicted by the DG approach 
are completely unphysical and much poorly agreeing with experiment that more common approaches, contrary to the 
seemingly original success claimed in Ref. |^. In Sec. VIII of Ref. [s], we clearly showed that a standard NEGF-DFT 
calculation yields currents slightly larger by a factor ^ 1.5 — 3, while DG's currents^ represent ~ 2 — 5% of the 
experimental currents of the recent accurate experiment of Ref. llOi 
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